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EFFECT OB' HIGH-VELOCITY FLUID FLOW OH THE 


BENDING VIBRATIONS AND STATIC DIVERGENCE 

OF A SIMPLY SUPPORTED PIPE 

By Harold I.. Dodds, Jr., and Harry L. Runyan 
Langley Research Center 


SUMMARY 




The present experimental study is concerned with the effect of high- 
velocity fluid flow cn the static and dynamic characteristics of a simply sup- 
port< d pipe. The frequency variation of the bending vibrations as a function 
of fluid velocity, as predicted by G. W. Housner, was verified experimentally. 
In addition, the critical fluid selocity at which the system becomes statically 
unstable was also verified. 
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INTRODUCTION 


High-perfornHu.ee liquid-prcr ellant launch vehicles require the rapid 
transfer of lurge quantities of fluid from tanks to pumps and engines through 
relatively thin-walled pipes. A study of the influence of the resulting high- 
velocity fluid on the static and dynamic characteristics of the pipes in launch 
vehicles is therefore necessary. 

Initial work in this field was motivated by vibration problems encountered 
in the transport of crude oil in pipes as well as possible problems in the feed 
lines of large water turbines. Ashley and Haviland (ref. 1) were \he first to 
investigate the bending vibrations of a simply supported pipe containing a 
flowing fluid. Hcusr.er (ref. 2) subsequently derived the equation of motion 
for the same case more completely and developed an equation relating the fun- 
damental bending frequency of a simply supported pipe to the velocity of the 
internal flow of the fluid. Housner also snowed that, at a certain critical 
velocity (similar to the Euler buckling formula for columns), a statically 
unstable condition could exist. Long (ref. 3) used an alternate solution to 
Housner' s equation of motion for tne simple end conditions and also treated the 
fixed-free end conditions. In addition, Long attempted to confirm the analysis 
by comparison with experiment. The experimental results were rather inconclu- 
sive since the maximum fluid velocity available for the test was low and the 
change in the bending frequency was very small. The maximum velocity obtained 
by Long was only 20 percent of the critical velocity which resulted in a 





maximum change in the bending frequency of only 3 • 2 percent. Other treatment-, 
of the subject were made by Ta Li (ref. 1), Niordson (ref. 3), and Benjamin 
(ref. 6). 

It is the purpose of this paper to report the results of an experiment to 
(a) study the effect of fluid velocity on the fundamental bending frequency of 
a simply supported pipe and (b) determine the existence of a statically 
unstable condition as predicted in reference 2. 

A number of engineering students, while participating ccnsecutivelv in 
the cooperative work study program at Langley Research Center, conducted the 
test program under the general supervision of the second author. The students 
who planned the experiment, selected and assembled the test apparatus, and 
conducted the tests were Messrs. Harold L. Dodds, Jr., Richard A. Crocker, 
Charles H. Fox, Jr., James E. Gardner, and A. David T.uckey. 
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amplitude coefficient in series expression (ea. (2)) for 

k = 1, 2, J, . . . 

amplitude coefficient in series expression (eq. (2)) for 
n = 1, 2, 3, . . . 

modulus of elasticity, lh/ft^ 

moment of inertia , ft 1 * 

integer 

length, ft 

mass of pipe plus fluid per unit length, slugs/ft 

integer 

time, sec 

average velocity of the fluid relative to x-direction, ft/sec 
coordinate measured along pipe 


coordinate measured perpendicular to pipe 

mass of fluid per unit length, slugs/ft 

angular velocity of pipe relative to its length, rad/sec 
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i mode of vibration, i = 1, 2, 5, . . . 

ANALYSIS 


The equation of motion of a uniform pipe containing flowing fluid is 
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which vac derived in reference 2 by using elementary beam c. cory. The coordi- 
nate system is shown in figure 1. The three terms on the right s; : j of 



Flour* 1.- Schematic d'jcram of exper' -ental setup. 


equation (l) represent the lateral inertial forces which act on the pipe. The 
first term is the force caused by the change in direction of velocity of the 
fluid due to the curvature of the pipe. The second term is the force associ- 
ated with the Coriolis acceleration which occurs because the pipe has an 
angular velocity with respect to its length while the fluid has a linear veloc 
ity with respect to the pipe. The last term represents the force which is 
related to the vertical acceleration of the pipe. For simplification, struc- 
firal damping has been omitted in equation (.1). 

For the case of free vibrations, the general solution of equation (l) for 
a simply supported pipe with normal modes z^ is given in reference 2 as 
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In order to evaluate the coefficients in this equation, equation (2) is sub- 
stituted into equation (1) and the resulting terns are grouped according to 
their trigonometric constituents. Equating the coefficients of like trigono- 
metric terms to zero in order to satisfy equation (l), the following set of 
equations for n = 1, 2, 3, . . . and k = 1, 2, 3, . . . are obtained: 
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Supressing all the modes except ch" first two, the simultaneous solution c:' 
equations ( 3 ) and (4) results in the following frequency equation for steady 
state vit rations: 
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Equation (5) gives the relation between V and for the pipe. 


By letting cu^ -» 0 for the case i - 1, equation (5) yields a critical 
velocity at which a static divergence occurs. 


V c 



(6) 


As pointed out in reference 2, replacement of pV^ in equation (6) by P 
provides an equation identical a the Euler buckling formula for columns where 
P represents the axial load. 


EXPERIMENTAL INVESTIGATION 


The apparatus and testing procedure which were used to study the effect of 
high-velocity fluid fl.w on the bending v .'orations and stability cf a simply 
supported pipe are described in the following section. 


Apparatus 

A schematic representation cf the experimental setup is shown in figure 1. 
The test pipe was constructed cf aluminum alloy (c0cl-T6) with a nominal out- 
side diameter of l.C-0 inch, a wall thickness of O.Co : j inch, and a length cf 
10 feet. The pipe w > 3 mounted horizontally with the downstream erd supported 
by two ball Leurin.v which allowed unrestrained horizontal translation and 
rotation of the d -.-.stream end. The upstream end was connected to a swivel 
joint. This type of support system provided an accurate simulation of a beam 
with pinried-pinned boundary conditions. The tall be-”-ir.gs at the end of the 
pipe prevented axial stresses during lateral deformations. A strain gage was 
mounted in the center of the pipe at the location of the maximum bending moment 
from which the frequency and damping of the bending vibrations could be 
determined . 

Tne average velocity of the fluid (water) was measured by means of an 
orifice plate located approximately 4 feet ups treats from the swivel joint. 

The instrumentation for measuring the fluid velocity consisted of the orifice 
plate, a differential-pressure transducer, and a galvanometer recorder. This 
equipment was accurate to within 2 percent (rms value) of full-scale response. 
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r The fluid for this experiment was obtained from «* hi gh -pres cure water 

I source located at the Langley landing loads track. This facility, capable of 

P delivering 750 gallons per minute at a maximum pressure of J200 lb/in2, is 

described in detail in reference 7* 


Testing Procedure 

In order to initiate the Tree vibrations of the pipe at various fluid 
velocities, a method of imparting a small disturb .nee to the pipe was necessary. 

This disturbance was accomplished by attaching a cord at the midpoint of the 
pipe which could be manually operated during tee test. The procedure for con- '$< 

ducting a test was to disturb the pipe by pin Ling the cord vertically, to 
record the response of the p spe, and to increase th' 1 fluid velocity continu- 
ously until instability of the pipe was achieved. 


DISCUSSION OF RESULTS 


This section includes both the calculated and experimental results and 
also seme remarks concerning the "damping" characteristics of the system. 

The values of the constants used in calculating the value of V c from 
equation (6) and the relation between and V frem equation (5) are as 

follows : 

p = 8.0C X 10”3 slugs/ft 
E = l.U x lo9 lbs/ ft 2 
I *• 1.0116 x 10“ 6 ft* 4, 
m = 1.19 x 10 -2 slugs/ft 
l - 10.5 ft 

A departure from an ideally uniform beam was required for the upstream 
end of the pipe in order to accommodate part of the swivel Joint (fig. 2). 

The Joint provided about 0.5 foot of rigid heavy pipe for an overall test 
length of 10.5 feet. The effect of this short stiff section ( compared with 
the thin-walled aluminum tube) on the critical velocity is considered negligi- 
ble since the curvature of the pipe is a maximum at its center and must 
approach zero at the ends in order to satisfy the boundary conditions. Cn the 
other hand, the effect of the swivel joint on the frequency of the pipe could 
he more pronounced because of the concentrated mass at the end of the pipe. 
However, comparison of the calculated frequency for a uniform pipe 10.5 feet 
long with the experimental frequency for the test pipe, which includes the 
deviation resulting from the swivel joint, indicates no appreciable effect on 
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Figure 2.- S.irtl |Oint at upstream end of test Qira. I 64-699L1 

the natural frequency. Consequently, the fundamental equations in refer- 
ence 2 were not rederived to include the possible effects caused by the 
swivel joints. 

The calculated results from equations ( 5 ) and (6) for the fundamental 
mode of vibration end the experimental data for two identical pipes are shown 
in figure 3- The specific experimental values pertaining to figure 3> which 
have been normalized by the theoretical values of cu at V = 0 and V c 
^a)y_o = 28. b rad/sec and V c = 129-5 ft/sec ), appear in the following table. 

The overall damping of the system at each dat3 point, as determined from the 
strain-gage records, also appears in the table: 



[ Log decrement | 


Not obtained Not obtained 
1*C53 I 0.5?9 
- 95 * -337 

.919 .629 

.676 .865 

.739 1.093 

0 Not obtained 


"*ror pipe 1 at v/v c = 0, the data vere 
obtained when the pipe va3 er*pty; however, 
the value, , t - ) v - o = 3. -0^2, hat been cor- 
rected to include the nast the fluid for 
a full pipe 
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Figure 3.- Comparison o! theoretical « id exoei imental resuits. 


As can be observed in figure 3> a good correlation between theory and 
experiment is achieved. The maximum experimental error reLative to the the- 
oretical curve is approximately 7<5 percent. 

As predicted by the theerj , the system became unstable which resulted in 
a permanent deformation of the pipe. Divergence was achieved with both pipes; 
however, because of an instrumentation failure during the first test-, only the 
divergence velocity for the second pipe was recorded which was 127-5 ft/sec. 
The experimental error of the divergence velocity relative to the predicted 
critical velocity is 1.5 percent. Pictures of the divergence phenomenon are 
shown in sequence in figure 4. Figure 4(a) shows the curvature of the pipe at 
V = 120 ft/sec which is just before the pipe diverged. Figure u(b) shows the 
curvature of the pipe at the divergence velocity of 127-5 ft/sec. Figure 4(c) 
shows the permanent deformation of the pipe alter the test was completed. 

The log decrement, as shown in the table, is a measure of the total 
damping of the system at each data point. These values were obtained by 
taking the natural logarithm of the ratio c f successive amplitudes for the 
first completed cycle of each strain-gage j -cord. The log decrement value 
when the pipe is empty can be considered a.v-a measure of the structural 
damping of the pipe and its end supports (ref. 5). This value was 0.179 for 
the first pipe. 

As may be observed in the table, the natuiai frequency of the pipe 
decreases as the fluid velocity increases wn lie the log decrement values 
increase. It should be noted that the log decrement increase is due to the 
combined effect of the frequency and fluid velocity, although the values dif- 
fer in some cases for the two pipes. 
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CCh'CIJTDXNC REMARKS 


The present investigation has been concerned with an experimental study 
of the effect of fluid velocity on the natural frequency of the free vibra- 
tions of a simply supported pipe. Tests were conducted witL two identical 
pipes both of which verified the predicted variation of the fundamental 
bending frequency with fluid velocity. The predicted critical fluid velocity 
for which the system becomes statically unstable was also verified experimen- 
tally. From the results presented, it appears that the theory of G. W. Housner 
is adequate for predicting the effect of fluid velocity on the natural fre- 
quency of simply supported pipes as well as the particular fluid velocity at 
which a static instability exists. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 15, 1565* 
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